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Abstract 

Effects of energy dissipation on quantum nucleation of two-ffavor quark matter 
in dense nuclear matter encountered in neutron star cores are examined at low tem- 
peratures. We find that low-energy excitations of nucleons and electrons reduce the 
nucleation rate exponentially via their collisions with the surface of a quark matter 
droplet. 
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The possible existence of quark matter in neutron stars has been studied for the past two 
decades from the viewpoint of equilibrium phase diagrams (see, e.g., Refs. 1)~3)). Recently, 
GlendenningB- re.axed the consent of local charge neutrality, as usuaHy aimed, and 
discovered that a phase where quark matter and nuclear matter in /3-equilibrium coexist in 
a uniform sea of electrons could intervene between bulk quark and nuclear matter phases 
for a finite range of pressures. This is because the presence of strange and down quarks 
plays a role in reducing the electron Fermi energy and in increasing the proton fraction of 
nuclear matter. GlendenningB and Heiselberg et al.& also claimed that the mixed phase 
could exhibit spatial structure such as quark matter droplets embedded in nuclear matter 
as a result of surface and Coulomb effects. 

As a star whose core consists of nuclear matter in /3-equilibrium spins down or accretes 
matter from a companion star, the central density could become sufficiently large for the 
mixed phase to be stable. Whether the mixed phase actually appears, however, depends 
on the occurrence of the dynamical processes leading to its appearance. The first of these 
processes should be deconfinement; the direct conversion of nuclear matter to three-flavor 
quark matter that is more stable than two-flavor quark matter is inhibited, since relatively 
slow weak processes are involved. If quark matter is deconfined from nuclear matter via a 
first-order phase transition, a droplet of two-flavor quark matter could form in a metastable 
phase of bulk nuclear matter during the spin-down or accretion. Iida and Sato!) calculated 
the rates for the formation of this droplet via quantum tunneling at zero temperature by 
incorporating the electrostatic energy into the potential for droplet formation. It was found 
that a droplet formed via quantum fluctuations would develop into bulk matter due to 
electron screening effects on the Coulomb energy of the droplet itself. This result implies 
that the mixed phase is unlikely to occur. 

The effects of energy dissipation on the quantum kinetics of the deconfinement transition, 
which were ignored in the previous work, B may have a decisive role in determining not 
only the degree of overpressure required to form an up and down quark matter droplet 
but also the crossover temperature from the nucleation via quantum tunneling to that via 
thermal activation; its comparison with the temperature of matter in neutron star cores is 
significant. In this paper, we thus estimate to what extent the rates for droplet formation 
via quantum tunneling are reduced by the friction exerted on the droplet, according to 
the work of Burmistrov and DubovskiiB that gave useful expressions for such reduction 
on the basis of a path-integral formalism developed by Caldeira and Leggett.S 1 It is found 
to be the ballistic motion of low-lying excitations of nucleons and electrons that controls 
the dissipative processes in the situation of interest here. The possible effect of neutron 
and proton superfluidity, which acts to decrease the number of available quasiparticle states 
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with momenta close to the Fermi surfaces, is also considered. Hereafter, we generally take 
h = c = ks = 1- 

The static properties of the deconfinement transition are described using simple phe- 
nomenological models for the energies of nuclear and quark matter, as in Ref. 4). For 
nucleons, we write the formula for the energy density as& 

e N (n,x) = n[m n + (K /18)(n/n - l) 2 + 5 (n/n ) 7 (l - 2x) 2 ] , (1) 

where m n is the neutron mass, n is the nucleon number density, n = 0.16 fm -3 is the nuclear 
saturation density, x is the proton fraction, K = 250 MeV is the nuclear incompressibility, 
and So = 30 MeV and 7 = 1 determine the symmetry term. For quarks we adopt an energy 
density based on the bag model, 

= (l |f) +B. (2) 

\q=u,d ) 

where q = u,d denote up and down quarks, n q is the number density of q quarks, and 
\i q = (1 — 2a s /n)~ 1 / 3 (7r 2 n q ) 1 / 3 is the chemical potential of q quarks. In Eq. (2), all quark 
masses have been taken to be zero, and the QCD fine structure constant and the bag constant 
have been set as a s = 0.4 and B = 120 MeV fm -3 . The energy density of the electrons, 
which are relativistically degenerate, is e e (n e ) = /Xg/47r 2 , with the electron number density 
n e and the electron chemical potential /x e = (37r 2 n e ) 1 / 3 . We neglect the energy contributions 
of muons and hyperons; the former was shown to be unimportant to the deconfinement 
transition itself, 0* whereas the latter will be considered elsewhere. Thermal effects on the 
static properties of the deconfinement transition may be safely omitted; for the temperature 
range considered, T < 10 MeV, we find T <C /Xj — rrii ~ 100-1000 MeV, where i = e, n 
(neutron), and p (proton), and /Xj is the chemical potential of i particles including their 
rest mass m^. The two electrically neutral bulk phases of interest are the /3-stable nuclear 
phase, in which /x„ = /x p + /x e , and the quark phase arising from the above-mentioned nuclear 
matter via deconfinement at fixed pressure P. The pressure Pq of the static deconfinement 
transition was obtained!- 1 from the comparison of Gibbs free energies per baryon between 
these two phases as Pq = 445 MeV fm -3 , where the baryon density jumps from ~ 5uq to 
~ 8n . 

The time r required to form a single quark matter droplet via quantum tunneling in the 
initial metastable phase consisting of nuclear matter in /5-equilibrium was evaluated & as a 
function of P within the WKB approximation & in the case in which no energy is dissipated. 
Let us now obtain a formula approximately reproducing the values of r so evaluated. Ac- 
cording to the results of Ref. 4), we expect that the contribution of the electrostatic energy 



Ec to the potential for droplet formation may be omitted by virtue of the electron screening 
effects on the droplet charge stemming from the pressure equilibrium between quark and 
nuclear matter. We thus write the Lagrangian for a fluctuation of the droplet radius R asB* 

L = (1/2)M(R)R 2 - U(R) . (3) 

Here M{R) is the effective mass derived from the nucleon kinetic energy as& 

M(R) = 4ire N {n init , x init )(l - n b , Q / n init ) 2 R 3 , (4) 

where n init and x init are the nucleon number density and proton fraction outside the droplet, 
and rib,Q = {n u + n d )/3 is the baryon density inside the droplet, and U(R) is the potential 
given by the minimum work needed to form the droplet asH* 

U(R) = (4irR 3 /3)n b , Q (p Q - // init ) + AE e (R) + A<naR 2 , (5) 

where /x ini t = [ejv(winit) ^imt)+e e (^e,init)+-P] / Minit is the chemical potential for the bulk nuclear 
matter with the electron number density n. e ,mit, — [^Q( n u,nd) + e e (n ei i n it) + P]/ n b,Q is 
the chemical potential for the bulk quark matter embedded in the electron gas of number 
density n e init , AE e (R) is the excess of electron energy over the initial value, and a is the 
surface tension which, being poorly known, is taken to be ~ 10 MeV fm -2 . Within the 
framework of the linear Thomas-Fermi approximation!^ in the determination of the electron 
distribution, which proves^ useful for R < 10 fm, we obtain AE e (R) = 47r/9Qi? 3 /i einit /3e, 
where pq = e[(2n u — rid)/ 3 — n et i nit ] is the excess charge density of the droplet, and /i e ,mit 
is the initial chemical potential of electrons. The quark number densities n u and rid are 
determined independently of R from the pressure equilibrium between quark and nuclear 
matter and from the flavor conservation before and after deconfmement. If one takes the 
energy Eq of the zeroth bound state around R = to be much lower than the barrier height 
Uq, an analytic formula for r is obtained as& 

t = ^ 1 exp(A ) , (6) 

where A is the doubled underbarrier action given by 

A = (hV2-K 2 /l6)\n biQ /nM t - M\f o-e N (n initl x ini t)R 7 J 2 (7) 

with the critical droplet radius R c = 3a/[nb t Q(pimt — Pq) — PQPe,mvt/ e ]i satisfying U(R C ) = 0, 
and u is the frequency of the oscillation around R = 0, Uq = [127r 3 / 2 r , (l/4)~ 2 ] 4 / 7 (47rcri? 2 ) 5/ ' 7 
x [M(R C )R 2 Y 2 I 7 . Naturally u is of order 10 23 sec -1 . As Fig. 1 illustrates, the values of 
lo 

Sio r gi ven by Eq. (6) at a — 10 MeV fm 2 agree fairly well with the computed values & in 
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Fig. 1. The formation time r of a single droplet with the surface tension a = 10 MeV fm 
as a function of pressure. The dashed line is the numerical WKB result 4 ) in the absence of 
dissipation; the dotted line is the result given by Eq. (6); the solid lines are the results for 
dissipative cases obtained from Eq. (9) as r = v^ 1 exp(A) at a = 0.2,1,5. Po denotes the 
pressure of the static deconfinement transition. 



the pressure range of interest. Discrepancies between these two results, which are shown to 
be ~ 20, arise mainly from the approximation E <C U = (A/21)AkgR 2 c made in Eq. (6) and 
partly from the approximation Eq = made in Eq. (5); the latter reduces the barrier height 
and width by < 10 %, and hence the order of the nucleation time by < 20. Hereafter, the 
Lagrangian (3) and the time (6) will thus be used as models for the formation of a droplet 
being in a completely coherent state. 

Complete coherence can be achieved only when all the excitations adjust adiabatically 
to the motion of the droplet surface. During the fluctuation of R, however, the low-energy 
excitations of nucleons and electrons appearing as quasiparticles in the vicinity of the Fermi 
surfaces undergo relatively slow relaxation processes, leading to dissipation of the total energy 
of the droplet. The fluid of quarks, which is assumed to be incompressible and hence static, ^ 
does not contribute to such dissipation. Since the mean free paths of i excitations (i = n,p, e) 
are estimated as k ~ n\ ^ 3 (uj — mi) 2 T~ 2 > 10 2 fm with the number density of i particles 
(n e = n p ~ no and n n ~ 4no), it is obvious that the k are large compared with the critical 
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droplet radius R c ranging 1-10 fm. Consequently, the excitations, which behave ballistically, 
collide with the droplet surface and the resulting transfer of momentum flux exerts an Ohmic 
frictional force on the droplet given byB ) 

F = -16tt - I)' {^a t n nm/lPF )j R 2 R , (8) 

where n nmji is the number density of a normal component of i particles, pp,i = (37r 2 ?T,j) 1//3 is 
the Fermi momentum of i particles, and is a factor of order unity that depends on the 
properties of i excitations and on their interactions with the droplet surface. For simplicity, 
we assume a e = a n = a p = a. It is instructive to note that F ~ —M(R)Rvp/R, where 

vf = (Ei "nm, iPF, 

i)/[e N (n init ,x init ) + e e (n e ,mit)] is the averaged Fermi velocity. 
We proceed to evaluate the effects of the dissipative processes mentioned above on the 
rates for droplet formation via quantum tunneling. First, we consider a system in which 
nucleons and electrons are in a normal fluid state at any temperature. The dissipation effects 
result in the multiplication of the tunneling probability exp(— Aq) by a reducing factor & 
exp(— AA); the expressions for A = Aq + AA were obtained!^ for an Ohmic frictional force 
such as Eq. (8) in the two limits A <C AA and A ^> AA. By noting that the typical 
velocity ^J2Uq/M(R c ) of a virtual droplet moving under the potential barrier is one order 
of magnitude smaller than vf (~ 0.4c), we find that the dissipation effects dominate the 
tunneling probability in the system considered here. The corresponding expression for A 
reads B* 

A = 4vra ( ^ - lV n am ,p F \ R 4 c s(T) , (9) 

V^init / \ i ) 

where s(T), the normalized underbarrier action composed of the potential term coming 
from Eq. (5) and the nonlocal dissipation term controlling the dynamics of the droplet in 
configuration space, is approximately 1.3 along the extremal trajectory for temperatures 
up to the quantum-classical crossover value T , as will be discussed below. The energy 
dissipation may also alter the oscillation frequency u by about one order of magnitude 
because of the difference between vf and \R\, but this has little effect. As depicted in 
Fig. 1, the nucleation time r = Uq 1 exp(A) derived from Eq. (9) decreases exponentially 
with increase in the overpressure AP = P — Pq, a feature resulting from the narrowed and 
lowered potential barrier. 

In a neutron star, the formation time of the first droplet is t/N, where iV is the number 
of virtual centers of droplet formation in the star. iV is roughly estimated to be V/(4ttRI/3) 
~ 10 48 , where V is the volume of the central region in which the value of r determined by 
the pressure remains within an order of magnitude, by calculating the pressure and density 
profiles of the star from the Tolman-Oppenheimer-Volkoff equations and the equation of 
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Fig. 2. The crossover temperature To from the quantum tunneling to the thermal activation regime 
as a function of a at a = 5, 10, 20 MeV fm~ 2 , N = 10 48 and t/N = 1 kyr. 



state given by Eq. (1). From these calculations we have confirmed that the central pressure 
of the star with maximum mass, ~ 730 MeV fm -3 , is high compared with the pressure range 
of interest. In order that the first droplet may form during a time scale of the spin-down or 
accretion T\ ~ yr-Gyr, the overpressure should amount to such a value AP C as to make t/N 
comparable to T\. Taking as an example a = 10 MeV fm -2 , a = 1 and N = 10 48 , we see that 
the dissipative processes change AP C from ~ 7 MeV fm~ 3 to ~ 21 MeV fm -3 and R c from 
~ 4.1 fm to ~ 1.9 fm. Correspondingly, the crossover temperature To from the quantum 
tunneling to the thermal activation regime, which is derived from the relation!) A = [To/To, 
goes from ~ 1.7 MeV to ~ 0.4 MeV. In Fig. 2 we have shown the crossover temperature T 
as a function of a at o = 5, 10, 20 MeV fm -2 and N = 10 48 . This figure suggests that the 
first-order deconfinement transition in nuclear matter of temperature < 0.1 MeV, as found 
in a neutron star core, if it occurs dynamically, is likely to proceed via quantum nucleation, 
although no definite conclusion can be drawn in the absence of the exact information about 
a and a. 

Finally, we consider the system in which neutrons (protons) are in a superfluid state for 
temperatures below the critical value T cn (T cp ). Whether or not T cn and T CJ , are finite 
for densities of about 5n is not well known, but it is expected from the properties of 
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nucleon-nucleon interactions that the 3 P2 neutron pairing and the 1 S'o proton pairing may 
possibly sustain a superfluid state (e.g., Ref. 10)). Such a state plays a role in reducing the 
number densities of the n and p excitations and thus in weakening the decelerating effects 
of the energy dissipation on the nucleation. In order to estimate the maximal influence of 
the superfluidity, we set n nmn = n nm , p = in Eq. (9), an assumption valid for T <C T cn 
and T <C T CjP . Accordingly, only the collisions of the electron excitations with the droplet 
surface contribute to the energy dissipation, leading to vf ~ 0.06c, which is comparable to 
the velocity of a droplet moving under the potential barrier. These considerations show that 
for a > 1 the nucleation time continues to be determined mainly by the dissipative processes 
rather than by the reversible droplet motion. For the foregoing parametric combination 
a = 10 MeV fm~ 2 , a = 1 and iV = 10 48 , we obtain AP C ~ 11 MeV fin" 3 , R c ~ 3.2 fm and 
T ~ 1.0 MeV. 

In summary, we have found that the effects of the energy dissipation on the formation 
of a u and d quark matter droplet via quantum tunneling in a metastable phase consisting 
of /3-stable nuclear matter, governed by the ballistic motion of the nucleons and electrons 
excited in the vicinity of the Fermi surfaces, act to increase the critical overpressure AP C 
considerably at low temperatures relevant to matter in neutron stars. It is to be noted that, 
for better estimates of this increase, density fluctuations associated with these excitations 
in the medium should be taken into account, since the hydrodynamic and thermodynamic 
descriptions of the kinetic and potential energies given by Eqs. (4) and (5) are taken for 
granted in the present theory of the quantum nucleation based on Ref. 9). The mechanism 
of the growth of the nucleated droplet in the pressure range P > P$ also remains to be 
examined; it may be controlled by the dissipative processes associated with the excitations 
in the hydrodynamic <C R) or ballistic (Zj ^> R) regime, or possibly by the detonation 
arising from the density jump due to deconfmement. Finally, we believe that the present 
analysis can be used for the study of the nucleation of quark matter in dense nuclear matter 
of temperatures ~ 10 MeV, as is of interest in stellar collapse. 
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